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Abstract 

Intensive search for a proper and realistic equations of state (EOS) is still continued for studying 
the phase diagram existing between quark gluon plasma (QGP) and hadron gas (HG) phases. 
Lattice calculations provide such EOS for the strongly interacting matter at finite temperature (T) 
and vanishing baryon chemical potential (//#). These calculations are of limited use at finite \ib 
due to the appearance of notorious sign problem. In the recent past, we had constructed a hybrid 
model description for the QGP as well as HG phases where we make use of a new excluded-volume 
model for HG and a thermodynamically-consistent quasiparticle model for the QGP phase and 
used them further to get QCD phase boundary and a critical point. Since then many lattice 
calculations have appeared showing various thermal and transport properties of QCD matter at 
finite T and ub = 0. We test our hybrid model by reproducing the entire data for strongly inter- 
acting matter and predict our results at finite ub so that they can be tested in future. Finally we 
demonstrate the utility of the model in fixing the precise location, the order of the phase transition 
and the nature of CP existing on the QCD phase diagram. We thus emphasize the suitability of 
the hybrid model as formulated here in providing a realistic EOS for the strongly interacting matter. 



PACS numbers: 12.38.Mh, 12.38.Gc, 25.75.Nq, 24.10. Pa 



I. INTRODUCTION 



Relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory (BNL), large 
hadron collider (LHC) at CERN and future compressed baryonic matter (CBM) experiments 
at GSI are providing a great opportunity to study the properties of strongly interacting mat- 
ter in the laboratory at extreme temperatures and baryon densities. The outcomes of these 
experiments will be of an enormous use in obtaining a proper equation of state (EOS) for 
high and low temperature QCD phases [1-3]. On the theoretical side, lattice calculations us- 
ing QCD thermodynamics provide a valid EOS for strongly interacting QCD matter but the 
method still lacks reliability for a matter possessing a finite density of baryons. Therefore, 
precise mapping of the entire QCD phase boundary and the location of a hypothesized QCD 
critical point (CP) still pose a challenging problem before the experimental and theoretical 
heavy- ion physicists today [4-6]. Lattice people use certain approximations in obtaining 
thermodynamical quantities at small baryon chemical potential and indicate the exis- 
tence of a cross-over chiral transition at fin = and T 170 MeV which ends at a CP for 
Ai c /3T C < 1.0 [7-8]. However, some recent lattice results showed doubts over the existence 
of a CP on the QCD phase boundary [9]. We still do not precisely know whether the con- 
jectured phase boundary is an outcome of deconfinement and/or chiral phase transition and 
whether CP appears only on the chiral phase boundary [10]. In such circumstances, it seems 
worthwhile to intensify the search for a realistic EOS which can suitably describe both the 
QCD phases. In this paper, our aim is to obtain such a description using phenomenological 
models for the QGP as well as HG phases separately and combine them suitably in a hybrid 
model formulation. We then proceed to test whether the model results exactly reproduce the 
most recent lattice data which have arrived for vanishing jj, B and at finite temperature and 
thereafter the precise location and the nature of the phase transition at CP are investigated 
in this model. 

Recently several attempts were made to explain the lattice data at low temperatures by 
using ideal hadron gas (IHG) description where interactions among hadrons are altogether 
neglected [11-15]. However, Andronic and collaborators indicated in a most recent paper that 
the features of the low energy lattice data arise if a hard-core repulsion among the various 
constituents of the hadron gas is incorporated as excluded volume effect [16]. Simultaneously, 
several papers addressed the high temperature lattice data by using a suitable quasiparticle 



description for QGP in which the constituents of QGP acquire a T and/or fx B - dependent 
mass [17-19]. These results suggest us to use a hybrid model type description in which low 
temperature phase can be described by a thermodynamically - consistent excluded volume 
HG model and high temperature QGP phase is suitably described by a thermodynamically 
consistent quasiparticle model. Recently we have used this hybrid model in constructing a 
first order deconfining phase boundary between HG and QGP by employing Gibbs' equilibria 
conditions and found that the boundary indeed terminates at a CP beyond which a crossover 
region exists [20, 21]. The main difference between our model and the one used by Andronic 
et al [16] is that we provide a hard-core size to baryons only and mesons in our model can 
overlap and penetrate into each other whereas Andronic et al have given the same hard 
core size to all the hadrons existing in HG. Recently we have demonstrated that our model 
successfully describes the multiplicity distributions, multiplicity ratios, shear viscosity to 
entropy density (r)/s) ratio, square of speed of sound (c 2 .) etc and our results compare well 
with the experimental data [22]. There are various thermodynamically consistent versions of 
quasiparticle model existing in the literature [17-19]. However, our quasiparticle description 
involves only two parameters i.e., A and To appearing in the effective coupling constant [21] 
in comparison to other approaches which have three or four parameters. 

We must emphasize the new results obtained by us in this paper in order to demonstrate 
its importance. In particular, we first calculate various thermodynamical quantities e.g., 
three times the normalized pressure density 3p/T 4 , normalized energy density e/T 4 , nor- 
malized entropy density s/T 3 , trace anomaly factor (e — 3p)/T 4 revealing a measure of the 
interaction present in the QCD matter and the normalized baryonic succeptibility (xf /T 2 ) 
etc for which the lattice data have recently become available at fi B = [11, 13, 17, 23-24]. 
We show their variations with temperature and compare our curves with the lattice data for 
the entire QCD matter. The comparisons thus demonstrate the validity of our hybrid model 
in providing a realistic description of both the phases of QCD matter. We also extend our 
studies of the above quantities for finite ji B so that our results can be put to test when the 
lattice calculations in future become feasible. Similarly we also show results for the ratio of 
pressure to energy (p/e), square of the speed of sound (cj.) and other transport ratios like 
shear viscosity to entropy density (rj/s) for the gluon plasma only and compare them with 
the available lattice results at /j, b — 0. Moreover, in this paper, we get the location of CP 
and determine the order of the phase transition by studying the difference in entropy density 



i.e., As/T 3 as well as difference in sound speed Ac 2 s between HG and QGP phases exactly 
at the phase boundary. We notice that these quantities vanish at the critical point and 
hence indicate a clear change in the order of the phase transition at CP. Recent studies of 
Csernai at al [25], Sasaki and Redlich [26] and Lacey et. al. [27] have revealed that rj/s ratio 
involves a cusp like feature in the graph near the critical point when we plot its variation 
with temperature T. Our results in hybrid model supports these findings and thus helps in 
precisely locating the critical point. 

The rest of the paper is organized as follows : In Sec. II, we give a brief outline for the EOS 
in QGP phase and a detailed prescription we have used to calculate various thermodynamical 
as well as transport properties. In Sec. Ill, we obtain the EOS for HG in the new excluded- 
volume model and derive the relations to calculate various thermodynamical and transport 
properties from it. In Sec. IV, we present a detailed comparison of our results with those 
from recent lattice calculations and also predict the variations of these quantities at finite 
[i B . Finally, Sec. V fixes the detailed comparisons of results regarding the order of the 
phase transition at CP on the deconfining phase boundary between HG and QGP and final 
conclusions are then presented. 

II. EOS FOR QGP 

The EOS for QGP in a quasiparticle framework as used in this paper has been described 
in detail in Ref. [28, 29] and the calculations regarding thermodynamical quantities like 
pressure, energy density, particle density etc. can be found in our earlier work [20]. Here we 
add how our prescription can be used to obtain the transport properties e.g., shear viscosity, 
speed of sound etc. In this model, we start with the definition of average energy density 
and average number density of particles and derive all other thermodynamical quantities 
from them in a consistent manner. The expressions for energy density, number density and 
pressure are [20]: 
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In Eq. (3), p(T, 0) is denned as follows [20]: 

£(^1 = |£ + /Vl(^). (4) 

In Eq. (1), €i(xil) = (x i l) 3 K 1 (x i l) + 3(xil) 2 K 2 (xil), where K x and K 2 are the modified Bessel 
functions with Xi — and index i runs for gluons, up-down quarks q, and strange quark 
s. Similarly in Eq. (2), Ii(xil) = (xil) 2 K 2 (xil). di are the degeneracies associated with the 
internal degrees of freedom and po is the pressure at T = T . Entropy density for QGP can 
be obtained using the expression: 

_ e + p - Ei Afrgj /rx 

where index i runs for different flavours of quarks. The baryonic succeptibility can be 
obtained from the following relation [30]: 

* 2 = ^ = diT q (6) 

Let us now calculate the transport properties. Our calculation for shear viscosity is based 
on the prescription used by Sasaki and Redlich [32] who have calculated shear as well as 
bulk viscosities for QGP in the quasiparticle model. By definition, shear (rj) and the bulk 
viscosities (£) are defined as coefficients of the space-space component of the deviations of 
the energy momentum tensor from equilibrium if a small perturbation causes the system 
to slightly deviate from its equilibrium value. Near equilibrium, one gets the total energy- 
momentum tensor as follows [31, 32]: 



T, u = T% + AT^, (7) 

here T® u = J dj0jz^-]f- [/ + /o] , is the energy momentum tensor of the system in equi- 
librium with /o (/o) as the equilibrium distribution functions for particles (anti-partciles) 
and d is the degeneracy factor of particles/anti-particles. Further, AT^ is the deviation 
in energy momentum tensor caused by a shift in the equilibrium. We can expand AT^ in 
time-time, space-time and space-space part as follows 



AT^ = AT 00 + AT 0i + AT i0 + A% (8) 



Now the space-space components of AT^ can be written as the sum of traceless and 
the scalar part [31, 32]: 

ATij = -C,5ijd m u n - rjWij, (9) 

which involve the transport coefficients eg, the bulk (() and shear(^) viscosities, respectively. 
In Eq. (9), W mn = (d m u n + d n u m — §5 m „<9iW l ) where u represents the flow velocity. Using 
relaxation time approximation, the shear viscosity in a medium composed of one type of 
particles/antiparticles can be obtained from the following expression [31] : 

1 f d 3 k k A 

V = ^fJ d T t/o(! ± /o) + /o(! ± h)] , (10) 

where k is the momentum and E = y/k 2 + M 2 is the energy with M being the thermal mass 
for quark or gluon, ± for fermion and boson (i.e., quark and gluon), respectively and fo (/o) 
stands for equilibrium distribution function for particle and anti-particle, respectively [31]: 

/o(/o) = (e ( ^ )/T ±l)- 1 , (11) 

where we use =p/x in f and / , respectively. In Eq. (11), we use ±1 for quark and gluon, 
respectively. In Eq. (10), r is the collision time. To calculate the collision time for quarks, 
antiquarks and gluons in QGP we use the following expressions [32,33]: 

T q{q) = T7\ ' ( 12 ) 
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Tg = Z-V^ • ( 13 ) 
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Here a s = g 2 /4:ir, is QCD running coupling constant where g 2 has been taken from Ref. [20] 
and N F is the number of effective flavour degrees of freedom. We get the final expression for 
shear viscosity by adding the contributions of all types of particles in Eq. (10). Similarly 
the speed of sound in QGP is : 

where p is the pressure and e is the energy density of the QGP as given in Eq. (3) and (1), 
respectively. 



III. EOS FOR HG 



Recently we proposed a new thermodynamically consistent, excluded- volume model for 
the hot and dense HG [20-22, 34]. Our approach incorporates the following new features. 
Besides thermodynamical consistency, our model uses full quantum statistics so that our 
model will remain valid even in the case of large //g. Moreover, we incorporated excluded- 
volume correction arising due to hard-core baryons. We further assumed that the mesons 
can overlap and fuse into one another and hence cannot generate any hard-core repulsion. 
This is one major difference between our model and other models [16]. Moreover, we have 
demonstrated [22] that our calculations give a very good fit to the experimental ratios of 
multiplicities and other transport coefficients. The grand canonical partition function in our 
excluded volume model of HG can be written as follows [20-22, 34]: 
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where is the degeneracy factor of ith species of baryons, E { is the energy of the particle 
(E{ = ^k 2 + mf), V t ° is the eigenvolume assigned to each baryon of ith species and hence 
Y^j NjV® becomes the total occupied volume where Nj represent the total number of baryons 
of jth species calculated in excluded-volume approach. 
We can clearly write Eq.(15) as: 

InZf* = F(l - V?)/^, (16) 
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where U represents the integral: 

_g^_ [°° k*dk 1 

and Aj = exp^) is the fugacity of the particle, nf is the number density after excluded- 
volume correction and can be obtained from Eq.(16) as : 



The total pressure of the HG in our model is [22, 35-36]: 
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In Eq. (19), the second term on the right hand side gives the total summed pressure from 
all the mesons here taken as pointlike particles. R = J2i n T^i gives the fractional occupied 
volume due to all types of baryons. The energy density of HG is obtained from the following 
relation : 

i ^ ' 3 

Similarly entropy density in our model is [22]: 

s = e HG + Phg ~ VBn B - ^ S n s ^ 

Similarly the baryonic succeptibility at /ib — can be calculated from the following equation 
[15] = 

XS = T 2 ( 9 V H T ) ■ (22) 
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Our calculation for the shear viscosity is completely based on the method outlined by Goren- 
stein et al. [37]. According to molecular kinetic theory, we can write the dependence of the 
shear viscosity as follows [37]: 

r/ocn/(|k|), (23) 

where n is the particle density, / is the mean free path, and hence the average thermal 
momentum of the baryons or antibaryons is: 

<' k » = V^A ' (24) 

and A is the Fermi-Dirac distribution function for baryons (anti-baryons). For the mixture 
of particle species with different masses and with the same hard-core radius r, the shear 
viscosity can be calculated by using equation [22, 38]: 
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where rii is the number density of the ith species of baryons (or anti-baryons) and n is the 
total baryon density. 



In order to calculate the speed of sound at constant s/n, we have used the method given 
by Cleymans and Worku [39]. The speed of sound at = is easy to calculate since it is 
sufficient to keep the temperature as constant [40, 41]. However, the speed of sound (c s ) at 
finite chemical potential can be obtained by using the extended expression 
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where the derivative dfis/dT and d\x s /dT can be evaluated by using two conditions, firstly 
by keeping s/n constant, and then imposing overall strangeness neutrality [22, 39]: 
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where L = nf + , is the sum of the strangeness density of baryons and mesons in the 
HG. Similarly K = nf + nf 1 , stands for the sum of anti-strangeness density of baryons 
and mesons. In all the above calculations, we have taken an equal eigen-volume V° = ^p^- 
for each baryon with a hard-core radius r = 0.8 fm. Many authors have used hard-core 
radii varying in the range of 0.5 fm and 1 fm [42]. However, Cleymans et al. have used 
r = 0.8 fm in their analysis [43]. It has also been found that proton has approximately 
exponentially decaying positive charge distribution with an effective mean radius of 0.8 fm 
[44]. Furthermore, if we consider the energy density existing inside a proton as 4£> where 
B is the Bag-constant, then we can write : AB = M p /V° where M p is the mass of a proton. 
Taking B 1 ^ = 170 MeV or B = 109MeV/fm 3 , we will get r w 0.8 fm. We have taken the 
contributions of all baryons and mesons and their resonances having masses upto 2GeV/ c 2 
in our calculation for the HG pressure. We have also used the condition of strangeness 
neutrality by putting ^ f Si(n* — nf) = 0, where Si is the strangeness quantum number of 
the ith hadron, and nf(nf) is the strange (anti-strange) hadron density, respectively [45]. 
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FIG. 1: Variation of 3 times of normalized pressure with respect to temperature at two different 
values of /ie = and 500 MeV in our hybrid model. Lattice data points at hb = are taken 
from Rcf. [11, 23]. Solid line shows our result in our excluded volume model of HG and dashed 
curve gives result in quasiparticle model at hb = 0. Dotted curve represents our HG result and 
dash-dotted curve shows our results in quasiparticle model, at hb = 500 MeV. Dash-tripple dotted 
curve presents the results obtained by Andronic et. al. by using RGSG model for HG [16]. 

IV. RESULTS AND COMPARISON WITH LATTICE QCD 

In Fig. 1, we have studied the variations of the quantity 3p/T A with temperature at 
two values of baryon chemical potential /ib = and /is = 500 MeV, respectively. At 
[is = 0, we have also shown the recent lattice results [11, 23] for comparison. We find that 
our results from the hybrid model are in excellent agreement with the lattice datapoints. 
For comparison, we have also given the results obtained by Andronic et al. [16] for the 
low temperature HG phase where they have used the excluded volume model of Rischke, 
Gorenstein, Stocker and Greiner (RGSG) and they have used a hard-core radius same for 
each hadron (i.e. meson and baryon) as r = 0.3 ± 0.05 fm. We notice that our curve shows 
much better agreement than their curve when compared with the lattice data. We also find 
that at hb — 0, the curve from our excluded- volume model is smoothly connected with 



14 
12 
10 
8 
6 
4 
2 




(e/T 4 ) 

LATTICE DATA 

(e/T 4 ) (u. =0) Our Model 

HG B 

(e/T 4 ) HG (n B =0) RGSG Model [16] 
(eAT 4 ) QGp (^ B =0) Our Model 

(e/T 4 ) (u. =500) Our Model 

HG b 

(eAT 4 ) QGp (n B =500) Our Model 



200 



250 300 

T (MeV) 



350 400 450 500 



FIG. 2: Variation of normalized energy density with respect to temperature at jjlb = and 500 
MeV in our hybrid model. Lattice data points at /j,b = are taken from Refs. [11, 23]. Shaded 
portion of a curve shows the calculation of Andronic et al [16]. Short-dashed curve with circular 
marker and long-dashed curve with circular marker presents our HG model results with hard-core 
radius (r) of each baryon as 0.6 fm at /j,b = and 500 MeV, respectively. 

the curve obtained in the quasiparticle calculation around T = 170 MeV. This result yields 
enormous faith in the use of the hybrid model and the values of the parameters in our model 
appear suitably adjusted. Here we stress that we are not using any additional parameters 
except those which were incorporated in our previous papers [20, 21]. 

In Fig. 2, we have plotted our results for the variations of normalized energy density 
(e/T 4 ) with respect to temperature at two values of \ib i.e., /ib — and /is = 500 MeV. For 
vanishing our hybrid model results again compare well with the lattice results. Here 
again the curves obtained in two models of the hybrid model are smoothly connected at 
T = 170 MeV but at finite \ib (= 500 MeV), we notice a discontinuity in the curves obtained 
for both the phases of QCD matter. This explains that at /j, b = MeV and T = 170 MeV, 
first-order phase transition does not occur whereas the discontinuity at /jl b — 500 MeV, 
indicates the presence of a latent heat in the transition. We have also shown the values 
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FIG. 3: Variation of normalized entropy density with respect to temperature at \x B = and 500 
MeV in our hybrid model. Lattice data points at fi B = MeV are taken from Ref. [11, 23]. Shaded 
portion of a curve represents the calculation of Andronic et al [16]. 

of e/T 4 at low temperature obtained by Andronic et al. by using RGSG model for HG 
phase and we observed that our result again shows better agreement with the lattice data in 
comparison to the results obtained by them. Furthermore, we have also shown the results 
obtained from excluded volume model for HG if baryons have a hard-core radius r = 0.6 
fm, in order to show the stability of our HG model results. We observe that the change in 
hard-core radius does not produce any noticeable change in the results. 

In Fig. 3, we show the variations of normalized entropy density (s/T 3 ) with temperature 
at two values of fi B (i.e., and 500 MeV), respectively. For fi B — MeV, our hybrid 
model yields good agreement with the lattice results. Here again the quantity s/T 3 is 
smoothly connected at vanishing \x B when we view the results of HG and QGP models but 
a disconnected graph appears at /i B — 500 MeV. Dashed-tripple dotted curve enclosing a 
shaded portion shows the results obtained by Andronic et. al [16] using an excluded volume 
model for the HG where the hard-core radius for each hadron (meson and baryon) is taken 
as r = 0.3 ± 0.05 fm and this gives rise to a shaded portion. At finite \x B (= 500 MeV), the 
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FIG. 4: Variations of trace anomaly (e — 3p) /T 4 with respect to temperature at /j,b = and 500 
MeV in our hybrid model. Lattice points at fiB = are taken from Refs. [13, 17, 24]. Shaded 
portion represents the calculation of Andronic et al. [16]. 

discontinuity in s/T 3 curves for both the phases shows the presence of latent heat involved 
in the phase transition from HG to QGP. 

In Fig. 4, we have plotted the results obtained for the trace anomaly factor (e — 3p) /T 4 
in our hybrid model calculations using HG and QGP equations of state separately at //# = 0. 
We further compare our results with the results obtained in a recent lattice calculation [13, 
17, 24]. We notice that our results yield an excellent fit to the lattice data. The success of 
our hybrid model which involves a separate and distinct description for both the phases (i.e., 
low temperature HG and large temperature QGP), is indeed excellent in reproducing the 
features of the lattice curves. We also show the results at finite baryon chemical potential i.e. 
He — 500 MeV and it shows that the peak of trace anomaly factor is found to shift towards 
the lower temperature side. One other important thing we observe is that by introducing 
the finite baryon chemical potential into the system the curves show a completely different 
behaviour at low as well as high temperature side. We also present here the result calculated 
by Andronic et. al. by dashed-triple dotted shaded curve [16] and we notice that our results 
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FIG. 5: Variation of normalized baryon number succeptibility with respect to temperature at 
Hb = in our hybrid model. Lattice data points at hb = are taken from Ref. [24]. 

yield better fit to the data. 

In Fig. 5, we show the variations of baryonic succeptibility normalized as xf /T 2 with 
temperature at \ib = and compare our results with the lattice QCD results. For //# = our 
hybrid model results again compare well with the lattice data points. Here again the curves, 
for x^/T 2 obtained for HG and QGP phases are smoothly connected at around T = 170 
MeV, which shows the presence of a cross-over transition between two phases. In Fig. 6, we 
show the variations of pressure to energy density ratio with respect to temperature at /ib — 
and we find that our results are well supported by the lattice QCD results. We can view a 
discontinuity in the behaviour of this curve when we go from low temperature phase to high 
temperature phase and the discontinuity and/or a minimum appears at around T = 170 
MeV. 

In Fig. 7, we show the variations of square of speed of sound (c 2 ) with respect to 
temperature at \ib = MeV and again a comparison is given with the lattice QCD results. 
For hb = 0, our results reproduces the lattice results very well. Here again the quantity 
c 2 shows a similar discontinuity and/or minimum where both the curves for both phases 
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FIG. 6: Variation oip/e with respect to temperature at hb = in our hybrid model. Lattice data 
points at hb = are taken from Refs. [11, 23]. 

i.e., HG and QGP join each other and, therefore it further supports a smooth crossover 
transition between two phases. We have also shown separately the curve for c 2 s obtained by 
Andronic et. al. by using RGSG model for HG phase only [16]. However, the features of 
the curve differ from the lattice data, although it also yields a minimum at around the same 
temperature. 



V. CRITICAL POINT (CP) AND ORDER OF PHASE TRANSITION 

In this section we attempt to show the use of hybrid model in constructing a deconfining 
phase boundary between HG and QGP by using Gibbs' criteria and fixing the precise location 
and nature of CP existing on this phase boundary. 

In Fig. 8, we show what will happen to a quantity depicting the change in the entropy 
density from the phase transition at CP of the phase diagram. We define the normalized 
difference j§ = (s/T 3 )qg P — (s/T 3 ) HG and demonstrate its variations with respect to the 
coordinates of the phase transition points lying at the deconfining phase boundary. We find 
that ^§ 0.0 and positive along the deconfining phase boundary in the case of first order 
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FIG. 7: Variation of speed of sound with the temperature at \ib = MeV as obtained in our hybrid 
model and compared with lattice results [11, 23]. Dash-tripple dotted curve shows the result of 
Andronic et al. [16]. 

phase transition which supports the absorption of a nonvanishing latent heat in the phase 
transition from HG to QGP. However, we surprisingly notice that j§ « exactly at the CP 
and thus CP can be taken as a point where the first order phase boundary terminates and 
phase transition changes its order. 

Sound velocity yields an important transport property of the QCD matter created in the 
nucleus-nucleus collision experiments because the hydrodynamic evolution of this matter 
strongly depends on it. Speed of sound is related to the speed of small perturbations 
produced in the QCD matter in its local rest frame. A minimum in c s has also been 
interpreted in terms of a phase transition point [39-40, 46-51] where a large number of 
degrees of freedom present in the medium causes difficulty in its propagation. Further, 
Chojnacki and Florkowski [52] have proposed that a shallow minimum in the speed of sound 
near the smooth joining of HG and QGP phases in a hybrid model description corresponds 
to the presence of a cross-over transition. They have obtained a temperature dependence 
in the sound velocity by exploiting QCD lattice simulation for the high temperature phase 




FIG. 8: Variation of (As/T 3 ) = (s/T 3 )qgp — (s/T 3 )hg with respect to coordinates of various 
phase transition points on the (T, he) phase boundary. We have used transition points from Ref. 
[20]. 

(i.e. T > 1.5T C ) and ideal hadron gas description at low temperature (i.e. T < 0.15T C ) and 
used further interpolations to connect them smoothly. In Fig. 9, we have separately shown 
the variation of square of speed of sound i.e., c 2 for HG and QGP media. We have also 
shown the difference {Ac 2 ) = (c 2 )qgp — (c 2 )hg and demonstrated its variations with respect 
to the coordinates of the phase transition points defining the deconfining phase boundary. 
We again find that Ac 2 ~ 0.0 at CP. Thus this result lends further support to our finding 
regarding the change of the order of the transition at the critical point. 

In Fig. 10 (a), we plot the variation of rj/s of gluon plasma with respect to T/T c , where 
T c is the value of the critical temperature. We compare our model result with the results 
obtained by Khvorostukin et. al. [53] and Bluhm et. al. [54]. A comparison with lattice 
calculations [55-56] is also demonstrated in Fig. 10 (a). We observe that the results obtained 
in our calculation agrees well with the lattice data even at very large temperatures. However, 
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FIG. 9: Variation of c^)qgp (dash-dotted line ), {c 2 )hg (dashed line), and (Ac^) = {c 2 s )qgp — 
((? s )hg (short-dashed line) with respect to coordinates of various phase transition points on the 
(T,/jlb) phase boundary. We have used transition points from Ref. [20]. 

the error quoted in the lattice simulation is quite large. Our result also appears in close 
agreement with the result obtained by Khvorostukin et. al [53] while the results of Bluhm et 
al [54] neither agrees with the lattice data nor shows any agreement with our calculations. 
It should be mentioned here that the authors [53-54] have used the QPM to calculate shear 
viscosity of a gluon plasma. The difference appears when they use the relaxation time (r) 
and its dependence on the strong coupling constant (g). Khvorostukin et al [53] have taken 
t oc g 2 while Bluhm et al [54] have assumed r oc G 4 where G is an effective strong coupling 
constant. This result naturally gives us additional confidence in the use of QPM for an EOS 
of QGP. In Fig. 10 (b), we plot the variations of shear viscosity to entropy density ratio 
for HG and QGP, separately with temperature at \Ib — 300 MeV (dashed-line) and at the 
critical potential fi c = 155 MeV (solid-line). At fig = 300 MeV, we observe a discontinuity 
in rj/s at the joining point of the curves for both the phases. Further we observe an upward 




FIG. 10: (a) Variations of shear viscosity to entropy density ratio (rj/s) for gluon plasma with 
respect to T/T c . Solid line is the result of our calculation. Long-dashed line is the result obtained 
in Ref. [53] and short-dashed line is taken from Ref. [54]. The lattice data with 16 3 x 8 and 24 3 x 8 
lattice are from Refs. [55] (triangles and squares) and [56] (filled circles), (b) Variations of i]/s 
obtained in our calculation for HG and QGP phases with temperature at different values of hb- 



jump in rj/s as we go from low temperature HG phase to high temperature QGP phase and 
this supports the result obtained by Sasaki and Redlich [26]. At critical chemical potential 
He = 155 MeV and temperature T c = 166 MeV [20-21], we get a cusp like behaviour in rj/s 
while going from HG to QGP phase as pointed out in Ref. [26, 57]. 

Thus above results give a firm indication that the order of phase transition changes at 
CP for the deconfinement phase transition. It should be added here that many authors in 
the past have also used two separate equations of state for QGP and HG and obtained a 
tentative explanation for an analytic and smooth cross-over and CP in their models [58-60] . 
Our model presents a similar picture. Here we explicitly and separately consider both the 
phases, i.e., HG as well as QGP and hence it gives a clear understanding how a first-order 
deconfining phase transition can be constructed in nature. At almost vanishing baryon 
density, overlapping mesons fuse into each other and form a large bag or cluster, whereas 



at high baryon density, hard-core repulsion among baryons restricts the mobility of baryons 
in the HG. Consequently we consider two distinct limiting regimes of HG, one is a meson- 
dominant while the other is a baryon dominant region and CP occurs exactly at the joining 
of two regimes in our models. 

Searching for the precise location and the nature of the critical point (CP) on the QCD 
phase diagram are still a challenging problem. Although various calculations have predicted 
its existence but the quantitative predictions regarding its location wildly differ. Experi- 
ments face an uphill task in probing the CP in QCD phase diagram because a clarity in 
theoretical prediction is missing. Moreover, many unstudied problems such as short lifetime 
and the reduced volume of the QGP formed at colliders also affect the location of CP and its 
verification [61]. In these circumstances, we hope that our results will clarify the mist sur- 
rounding the understanding of the deconfining phase transition. More importantly, we have 
formulated a phenomenological hybrid model which provides a realistic EOS for the entire 
QCD matter and in the absence of the first-principle lattice QCD calculation especially 
at finite it can be reliably used for deriving the information on the QCD phase boundary. 
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